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Abstract 

In a companion paper (|hep-th/0512317 1, we have presented an approximation scheme to solve 



the Non Perturbative Renormahzation Group equations that allows the calculation of the n-point 
functions for arbitrary values of the external momenta. The method was applied in its leading 
order to the calculation of the self-energy of the O(A^) model in the critical regime. The purpose 
of the present paper is to extend this study to the next-to-leading order of the approximation 
scheme. This involves the calculation of the 4-point function at leading order, where new features 
arise, related to the occurrence of exceptional configurations of momenta in the flow equations. 
These require a special treatment, inviting us to improve the straightforward iteration scheme 
that we originally proposed. The final result for the self-energy at next-to-leading order exhibits 
a remarkable improvement as compared to the leading order calculation. This is demonstrated by 
the calculation of the shift ATc, caused by weak interactions, in the temperature of Bose-Einstein 
condensation. This quantity depends on the self-energy at all momentum scales and can be used 
as a benchmark of the approximation. The improved next-to-leading order calculation of the self- 
energy presented in this paper leads to excellent agreement with lattice data and is within 4% of 
the exact large N result. 
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I. INTRODUCTION 



The development of non-perturbative methods is essential to be able to deal with a 
large variety of problems in which the absence of a small parameter prevents one to build 
solutions in terms of a systematic expansion. Among such methods, the non perturbative 
renormalization group (NPRG) [3,0,0,0,0 stands out as a very promising tool, suggesting 
new approximation schemes which are not easily formulated in other, more conventional, 
approaches in field theory or many body physics. The NPRG has been applied successfully to 
a variety ofphysical problems, in condensed matter, particle or nuclear physics (for reviews, 
see e.g. 0,0, 0])- In most of these problems however, the focus is on long wavelength modes 
and the solution of the NPRG equations involves generally a derivative expansion which only 
allows for the determination of the n-point functions and their derivatives at small external 
momenta (vanishing momenta in the case of critical phenomena). In many situations, this 
is not enough: a full knowledge of the momentum dependence of the correlation functions 
is needed to calculate the quantities of physical interest. 

For this reason, in ref. we have presented an approximation scheme to solve the 
NPRG equations that allows the calculation of the n-point functions for arbitrary values of 
the external momenta. The method was applied in its leading order to the calculation of 
the self-energy of the O(A^) model in the critical regime. The purpose of the present paper 
is to extend this study to the next-to-leading order of the approximation scheme. This 
involves the calculation of the 4-point function at leading order, where new features arise. 
In particular we need to deal with exceptional configurations of momenta that enter the flow 
equations. Because of these, the straightforward iteration scheme proposed in ref. Q yields 
some unphysical features in the 4-point function. After having identified the origin of the 
problem, we shall show how it can be cured by a proper treatment of the flow equation in 
the channel where the exceptional momenta matter. The final result for the self-energy at 
next-to-leading order exhibits a significant improvement as compared to the leading order 
calculation. In particular the calculation of the the shift ATc in the transition temperature of 
the weakly repulsive Bose gas a quantity which is very sensitive to all momentum scales 
and which is used as a benchmark of the approximation, is now in excellent agreement with 
the available lattice data, and within 4% of the exact large N result. Note that preliminary 
results concerning the calculation of ATc at next-to-leading order have been presented in 
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ref. These results were obtained without the improvements just alluded to. 

This paper is a sequel of ref. Q], and should be read in conjunction with it (hereafter 
ref. will be referred to as paper I, and the prefix I in equation labels will refer to equations 
in paper I). As in paper I we shall focus the discussion on the O(A^) model, although most 
of the arguments have a wider range of applicability. Thus, we shall consider a scalar ip^ 
theory in d dimension with 0{N) symmetry: 

^ = J^'^{1 + ^ + ^ , (1) 

The field ip{x) has real components ipi{x), with i = 1, - ■ ■ ,N. 

The basic equation of the NPRG is the flow equation for the effective action Tf^[(j)] which 
interpolates between the classical action S and the full effective action r[0] (0 is the expec- 
tation value of the field) as the parameter k, varies from the microscopic scale A down to 
zero. The flow equation for the effective action Tf^lcp] reads 0,0,0,0,121: 

d.TM = / (1^ [rf + R.] , (2) 

where the trace runs over the O(A^) indices, T^^^ is the second derivative of with respect 
to 0, and Rniq) is the regulator chosen, as in paper I, of the form [l2| 



RM^i'^'-qyi^'-Q')- (3) 

The role of the regulator is to suppress the fluctuations with momenta q ^ k, while leaving 
unaffected those with q ^ k. 

The flow equations for the various n-point functions are obtained by taking functional 
derivatives with respect to of the equation for T^lcj)] (see eq. (1.8)). In particular, the 
self-energy S(/s:;p) is obtained by integrating the flow equation (1.9): 

dn^fh^^'^P) = Si2d^J:{K;p) = "2 y J^:^9^Ri^{q)G^{K;q)r'f]ii{K;p,-p,q,-q), (4) 
with the inverse propagator given by 

G-\K,q)=q' + E{K;q) + RM- (5) 

In eq. (^, and later in this paper, we often denote the O(A^) indices simply by numbers 1, 2, 
etc., instead of ii,i2, etc., in order to alleviate the notation. 
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The flow equations for tlie n-point functions constitute an infinite liierarcliy of coupled 
equations (for example, eq. (j3)) for the 2-point function contains in its r.h.s. the 4-point 
function). In paper I we have proposed a strategy to solve this hierarchy by following an 
iterative procedure. This starts with an initial ansatz for the ra-point functions to be used 
in the right hand side of the flow equations. Integrating the flow equation of a given n-point 
function gives then the leading order (LO) estimate for that n-point function. Inserting 
the leading order of the n-point functions thus obtained in the right hand side of the flow 
equations and integrating gives then the next-to-leading order (NLO) estimate of the n-point 
functions. And so on. 

Recall that there is no small parameter controlling the convergence of the process, and the 
terminology LO, NLO, refers merely to the number of iterations involved in the calculation 
of the n-point function considered. Since the calculations become increasingly complicated 
as the number of iterations increases, the success of the procedure relies crucially on the 
quality of the initial ansatz. A major task then is to construct such a good initial ansatz. 

The equations are solved starting at the bottom of the hierarchy, that is, with the equation 
for the 2-point function which involves, in its right hand side, the 2-point function (through 
the propagator), and the 4-point function. As initial ansatz for the propagator, we take the 
propagator of a modified version of the derivative expansion that we called the LPA'; it is 
given by (see paper I, sect. II): 

GllA'i^^'^fl) = Z^q^ + ml + R^{q), (6) 

where the field renormalization factor and the running mass m^^ are obtained by solving 
the LPA' equations (see paper I, sect. IIB). The initial ansatz for the 4-point function 
is given explicitly in paper I, sect. Ill, and is obtained as the solution of an approximate 
equation (see eq. (1.68) and eq. below). For more clarity, we shall distinguish in this 
paper the initial ansatz for F*^'*) by a tilde. For consistency, we shall use a similar notation 
for the initial ansatz for the propagator, i.e., we shall set G = Glpa'- Summarizing, the 
leading order self-energy is given by: 

5i2<9«Sio(ft:;p) = j -^^d^R^{q)G'^{K]q)f%{K]p,-p,q,-q). (7) 

Note that, as compared to eq. (jll), eq. ((Zj) is now a trivial flow equation since all quantities 
in the r.h.s. are known quantities: Sio('t;p) is simply obtained by integrating the r.h.s. 
with respect to k. The leading order self-energy has been studied in detail in paper I. 
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As stated earlier, the purpose of the present paper is to calculate T,nlo{p), the self-energy 
at next-to-leading order. To do so, we need to use in the r.h.s of eq. (j3)) the leading order 
expressions for both the propagator and the 4-point function. The leading order propagator 
is obtained from eq. (0) with the self-energy S^o- Getting the leading order expression for 
the 4-point function T^^^h{k;p, —p, q, —q) will be the main task of this paper; it is presented 
in sect. [H] First, in sect. Ill Al we follow the procedure outlined above, i.e., replace in the 
r.h.s. of the flow equation for the 4-point function, eq. (jH)) below, the initial ansatz for the 
propagator, the 4- and the 6-point functions. The initial ansatz for the 6-point function is 
obtained by following the same strategy as that used in paper I, sec. Ill, in order to construct 
the initial ansatz for the 4-point function. Although conceptually straightforward, this is 
technically more involved and the details are presented in app. EI Then, in sec. Ill Bl we 
present an improved procedure to calculate Tf2ii{K\ p, —p, q, —q) at LO, which copes properly 
with the difficulties related to the exceptional configuration of momenta. The properties of 
Satlo are presented in sect. IIIH together with the result of the calculation of A((y9^), which 
we use, as we have recalled earlier, as a benchmark of the approximation scheme. The last 
section summarizes the results, and points to further improvements of the approximation 
scheme that we have already started to implement HQ. 



II. THE 4-POINT FUNCTION AT LO 



The fiow equation for the the 4-point function in vanishing field reads (see e.g. eq. (1. 11)): 

dKTY2ii{fi]P, -P, g, -q) = J j^-^d^R^{q')G'^{t^] q') 

X [g{k] q')Tfl.{K] p, -p, q', -q')T[fl{K] q, -g, -q', q') 

q' + P + q)Tflj{n; p, q, q' , -q' - p - q)T2ilj{i^; -p, -q, -q', q' + P + q) 
+G{k; q' + p- g)rSJ^(K; p, -q, q', -q - p + q)T^r2ij{'^'^ 9' ~P^ "^'^ - ? + P) } 

We have specified here the O(A^) indices and the momenta of the particular 4-point function 
v'f^iiii^] P, —p, q, —q) that is needed for the calculation oiY^NLoip) (see e.g. eq. (jU). Following 
the terminology introduced in paper I, we refer to the second line of eq. (jH)) as to the s- 
channel contribution, while the third and fourth lines are respectively the t and w-channel 
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FIG. 1: Diagrammatic illustration of the contribution of the s-channel to the flow of the 4-point 
function (the second line in eq. ©)• The dotted line indicates the loop integral involving this 
contribution of F^^^ in the calculation of the self-energy in NLO, i.e., T,'^^^q{p) (see sect. III). 




FIG. 2: Diagrammatic illustration of the contribution of the t and u channels to the flow of the 
4-point function (third and fourth lines in eq. ((SI)). The dashed line indicates the loop integral 
involving this contribution of F^^^ in the calculation of the self-energy in NLO, i.e., sj^^Q(p) (see 
sect. III). 

contributions. Note that the contribution of the w-channel differs from that of the t-channel 
solely by the change of sign of q. A graphical illustration of these various contributions is 
given in figs. [T]|ni 

As stated above, the leading order expression for F*^^^ is obtained by substituting in the 
r.h.s. of eq. (jH)) the initial ansatz for G, F^^^ and F'^^^ In fact, we shall proceed with a further 
simplification which consists in setting g' = in the vertices of the r.h.s. of eq. (jH)). For F^^-*, 
this is justified by the fact that the initial anstaz F*^^^ varies little when the momenta on 
which it depends are varied by an amount smaller than k. This property has been explicitly 
assumed in the construction of F*^^^ in paper I. It has also been tested quantitatively in the 
calculation of the leading order self-energy (see fig. 15 in paper I and the discussion at the 
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FIG. 3: Diagrammatic illustration of the contribution of the 6-point function to the flow of the 
4-point function (last line in eq. ©)• The dotted line indicates the loop integral involving this 
contribution of F^^^ in the calculation of the self-energy in NLO, i.e., S^^^q(p) (see sect. III). 

end of sect. IV of paper I). We assume that this property also holds for the initial ansatz 
i'uiimmi'^'-'P^ — p, q, —q, q', — g')- ia.ct, for this latter function we shall also set g = 0, which 
can be justified as follows. Note first that T^^\p, —p,q, —q) will eventually be used in the 
calculation of T,nlo{p)i and in this calculation g ^ k. As we explained in paper I, sect. IIIA, 
setting g = is then well justified when p ^ k, because in that case all the momenta are 
smaller than n, and F*^^-' is well approximated by the LPA'; it is also justified when p ^ k 
since then q is negligible compared to p. It is only in the small integration region q k 
that the approximation could be less accurate. Observe finally that the contribution of F*^^) 
is negligible unless p <^ Kc- Thus, in line with approximation Ai of paper I, we shall, in the 
r.h.s. of eq. (jH)), set g' = in the 4-point functions F*^'*-' and g = g' = in F*^^). We then 
arrive at the following simplified equation: 

^^^^rSg^/t; p, -p, g, -g) = lf\K)r''^i-{K; p, -p, 0, 0)fJfA(fi;; g, -g, 0, 0) 
+jf\K;p + g)f S;^(k;p, g, 0, -p - g)f^Sj (/t; -p, -g, 0,p + g) 
+jf\K;p- q)t^^l{K;p, -g, 0, -p + g)f ^|^(k; -p, g, 0,p - g) 
-1/f («:)f!?,_(«:;p,-p,0,0,0,0), (9) 

where the functions I^J^\k) and J^^\k]p) are defined in eqs. (1.42) and (T55), respectively. 
The construction of the initial ansatz ff2iimmi^'^P^~P^^^^^^^^) requires the solution of an 
approximate flow equation which is obtained by following the same three approximations 
that are used in paper I, sect. Ill, to get r*^^^ This is presented in app. ^ The explicit 
traces of products of the functions F^^^ appearing in eq. ^ are given in app. B. 
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Eq. (0) will be used to calculate T^^^ at LO. Note that, as was the case for eq. ((Tj) giving 
Slo; eq. is now a trivial flow equation: all quantities in its r.h.s. are known quantities. 
We shall refer to this calculation of T^'^^ at LO, which follows strictly the scheme proposed in 
paper I, as to the "direct procedure". The results obtained in this way are discussed in the 
next subsection. We shall see there that this procedure yields unphysical features in some 
specific situations, whose origin will be discussed. An improvement on the direct procedure 
will then be proposed in the following subsection. 

Before proceeding further, let us mention that we have used the LO estimate of 
-'^i2h('^'P' ~P^^^ ~i) obtained in the direct procedure to perform a consistency check of ap- 
proximation Ai used in order to obtain eq. Q. We have verified that T^^^ui'^yP^ ~PjQj ~q) 
varies little as q varies in the range q < n, which is the range relevant for the calculation of 
^NLo{p)- Only in a small region where k is of order p <C Kc, can the function change by as 
much as 10% when q goes from to k. In all the other regions the variation is less than 1%. 
In the calculation of A(y9^) that will be reported in the next section, one needs J^Nioip) for 
values of p around k^. there, the approximation Ai is indeed excellent. However, at very 
small momenta, the error due to this approximation on the magnitude of T^j^iioip) can be 
large. But in this region the approximation Ai is not the dominant source of error anyway: 
in particular, any error on the exponent rj will translate into a large relative error on the 
magnitude of T,nlo{p)- 

A. Direct procedure 

As we have just discussed, since all quantities in the r.h.s. of eq. are known, 
T^^2iii'^'iP^ — p, g, —q) at LO is obtained by simply integrating the r.h.s. of eq. between A 
and K, and adding the bare value of T^^^ (i.e., the value of F*^^^ at the microscopic scale A): 

r%{^ = A; p, -p, q, -q) = {N + 2) 6,2 (10) 

with g\ = u/3, u being the parameter of the classical action (P) (see paper I, sect. IIB). 

The LO value of T^^2iii'^^P^ ~P^ O5 0) thus obtained is compared with both the initial ansatz 
and the LPA' result, (A^ + 2)gt^, in fig. 0] below. When k, ^ p one expects the LPA' to be 
a good approximation, and indeed the three curves almost coincide (the initial ansatz is by 
construction identical to the LPA' for large k) . When k goes to values smaller than p, the 
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FIG. 4: The 4-point function T^^\k;p, —p, 0,0) (in units of A) as a function of k/u: the four 
curves represent respectively the LPA' (full line), the initial ansatz (dotted line), the direct LO 
(dashed line) and the improved LO (dot-dashed line). The calculation is done for (a): p = u, and 
(b): p = u/10. 

flow of the 4-point function is expected to be slower than that of the LPA' since, generally, 
momenta in the propagators tend to suppress the flow. Fig. lU shows that it is indeed the 
case (in the initial ansatz this feature is implemented in a sharp manner with the help of 
6-functions depending on a parameter a (see paper I, sect. Ill, and eq. (fTTj) below)). As 
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shown by fig. El as n/p decreases further, the LO curve remains close to that corresponding 
to the initial ansatz but when n/p becomes too small, eventually the two curves deviate: 
while the initial ansatz goes to zero as k — 0, the LO curve goes to a negative value. As we 
shall explain shortly this is an unphysical feature of the LO in the direct procedure. 

At this point, it is instructive to recall the form of the approximate fiow equation whose 
solution is the initial ansatz V^^\ This equation is established in paper I, sect. Ill, and reads 
(see eq. (1.68)): 

x{fSti(P, -P,0,0)fg(g,-g,0,0) 

+ Q{k^~ a\p + qf) V%{p, q, 0, -p - q)T%i~p, -g, 0,p + q) 
+ Q{k'- a\p - qf) V%{p, -g, 0, -p + g)fg^(-p, g, 0, p - g) } . 

(11) 

The function is defined in eq. (1.44). It measures (approximately) the magnitude of 
the contribution of the 6-point function relative to that of the terms containing the 4-point 
functions. Note the similarity between the r.h.s. of eq. (jH)) giving the LO expression of r(4) 
and the r.h.s. of eq. (jlip for the initial ansatz V^^^: the main differences are the replacement 
of the O-functions of eq. (fTT|) by the proper loop integrals in eq. (jHl), and a more accurate 
treatment of the contribution of the 6-point function in eq. (jH)). This observation leads 
us to expect that F^^^ at LO should not differ much from T'^'^\ a necessary condition for 
the validity of the iteration procedure. As shown by fig. 01 this expectation is fulfilled, 
except for small values of k: while the LO F^^^ goes to a negative value when k — 0, the 
initial ansatz V'^^\k]p, — 0, 0) vanishes (one reads from eq. (1.99) that in the limit k ^ 0, 
f *^^^(k; — p, 0, 0) goes to zero like (^i^^^''^*'^^^''). There is in fact a major difference between 
eq. flllj) for the initial ansatz and eq. giving F^^^ at leading order: while in eq. (fTT|) the 
unknown function r*^^''(K; p, — 0, 0) sits in the r.h.s, this is not so in eq. Q, as we have 
already emphasized. Thus the structures of eqs. (fTT|) and Q are different, and only eq. (fTT|) 
captures the essential feature that guarantees the vanishing of V'^^\k]p^ —p, 0, 0) as k — > 0, 
which, as we shall show now, is a property of the exact solution of the flow equation. 

To do so, we shall consider eq. (jHl) for T'^'^^n^p, — p, 0,0) and study the behaviour of its 
solution in the limit when k 0. Observe flrst that the s-channel controls the flow when 
K ^ p: indeed, in this channel, the exceptional conflguration of momenta (p, —p) makes 
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the loop integral in the r.h.s. of the flow equation independent of p. This is manifest in 
eq. 0: the loop integral reduces to the momentum independent function while 
in the other channels the p-dependent functions J^'^^ {p ± q) enter (this is also obvious in 
eq. (jllj) where the functions J^\p ± q) are approximated using G-functions). It follows 
that in the s-channel, the external momentum p does not contribute to stabilize the flow 
whenever it does in the other channels: the flow continues all the way down to 

K = 0. When treated correctly, this is what induces the vanishing of T^^\k]p, — p, 0,0). To 
show this latter point, we focus on the contribution of the s-channel and, in the r.h.s. of 
eq. we keep T^'^\k;p, —p, 0,0) as an unknown function (instead of replacing it by the 
initial ansatz T^^\k;p, — p, 0,0)). We also temporarily neglect the contribution of the other 
two channels, and also the contribution of the 6-point function. Then, eq. (jH)) takes the 
following simple form (here we focus on the structure of the equation, dropping the O(A^) 
indices for simplicity, as well as explicit reference to the external momenta; a more precise 
equation is written in the next subsection): 

Kd^r{t,) = gji'\t,)r{^), (12) 

where we have replaced r*^'^^(K; 0, 0, 0, 0) by the LPA vertex g,^ (see eq. fll9p below). In the 
limit K, ^ (see paper I), (yf^ ~ Z"^ k,'^~'^ and ~ k'^~~'^/Z'^, so that g^lf\K) C,, where 

^ is a positive constant. It then follows form eq. ()12|1 that ~ as k ^ 0, and hence 
vanishes as n vanishes. Turning now to the contributions of the t and -u-channels, we note 
that these are more regular than the s-channel contribution (this can be seen for instance 
from the fact that the ratio jf\K]p)/lf\K) ^ as k/p ^ 0; see paper I, fig. 9). Similarly, 
the contribution of the 6-point function is proportional to /^^^ and / ijj^^ ~ k^"^'' as 
K — > 0. Using these properties, one can write the equation for the 4-point function at small 
K in the schematic form 

Kd^T{K)=^T{K) + ^K), (13) 

where $(k) may be considered at this point as a known function of n ($(k) is constructed 
from the initial ansatz for the 4-point and 6-point functions). The equation above can be 
easily solved, with the result 



T(k) 
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Since, as we have just argued, vanishes as k ^ 0, the integral does not diverge faster 
than K~^, and the small k behavior is governed by the factor outside the brackets, that is 
by the solution of eq. (jl2p . which guarantees that r(0) = 0. This argument shows also that 
provided one solves consistently eq. (jl3|) small errors in the estimate of the function $ are 
damped by the factor (k/kq)^. In particular, the solution for the initial ansatz is identical to 
that written above, eq. (fT^ . with $ replaced by another function $ not too different from 
$ (we may ignore here the factor F{k,) in eq. (fTT|) . which complicates the analysis, but in an 
inessential way. Accordingly we can also ignore the contribution of the 6-point function in 
eq. 0. Then the only difference between eqs. (jllll and (jHI) comes from the approximation 
of the functions J^^^ of eq.Q with the functions 9 in eq. (fTT|) ). One then expects the two 
solutions corresponding to $ and $ to be close to each other. 

Let us however imagine that we apply the direct procedure to our schematic equations, 
by caculating F^o from the analog of eq. Q, namely 



where the integral is convergent and yields a finite value for Tlo{i^ = 0). Thus, although the 
two solutions T{k,) and T{k,) would differ little if $ and $ differ little, the direct procedure 
which consists in simply integrating the r.h.s. of the flow equation fails to reproduce the 
low K, behavior. This small k behavior can only be obtained if the feedback of the flow is 
properly taken in the solution of the flow equation. But, as revealed in the previous analysis, 
this needs to be done only in the channel where the flow is not controlled by the external 
momenta, i.e., the s-channel. We shall implement this improved strategy for the LO in the 
next subsection. 

B. Improved approximation for r^^i^j^(K;p, —p,0,0) 

As we have seen, the main problem with the direct procedure that we have followed in 
the previous subsection comes from the s-channel where the exceptional configuration of 
momenta does not contribute to stabilize the flow when k — > 0. In this subsection and 



(15) 



where we have used the fact that K,di^r{K,) = ^r(/t) + $(/t). We would then obtain 




(16) 
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the following one, we shall present a more accurate treatment of this particular situation, 
making more precise the treatment presented in the previous subsection. 
To this aim, we consider first the case q = and replace eq. (jH)) by: 

Kd,r^S,,{K- p, -p, 0, 0) = Jf (/t)r;S,(/t; P, -P, 0, 0)f 0, 0, 0, 0) + $i2(/t; P), (17) 

where we have isolated the contribution of the s-channel, and grouped the other contributions 
into the function $i2(k;p): 

^i2if^,p) = jPiK;p)t%iK;p,0,0,-p)t%iK;-p,0,0,p) 
+ /f\K;p)t%{K;p, 0, 0, -p)T%{>^^ -p, 0, 0,p) 

- li!i\^)^fiimJ^;p,-p,W,o). (18) 

$i2(/«,p) is a known function which involves the initial ansatz F*^^^ and T^^\ In contrast, 
r[2;;(/t;p, — p, 0, 0) is taken to be the same function in the r.h.s and the l.h.s. of eq. (fTTjl . 
That is, eq. ()17|) properly takes into account the feedback of the flow in the s-channel on 
the solution of the flow equation. Note that for p = the solution of eq. (fT7|) (or eq. ©) is 
the LPA solution. We can therefore replace in the r.h.s. of eq. f|T7jl f j^^].;(/t; 0, 0, 0, 0) by the 
LPA value 

f[%iiK; 0, 0, 0, 0) = giK){5,,6M + 6,k6ji + 6u6,k)- (19) 

Then eq. (fTTj) becomes 

Kd^T^^\K;p) = {N + 2)I^'\K)g{K)T^'^\K;p) + ^k;p), (20) 

where we used the definitions 

rS,(«:;p,-p,0,0) = 5i2r(^)(/s:;p), <^i2{^;p) = 5i2$(«:;p), (21) 

and we omitted a factor 612 in both sides of eq. (PUI) . 

Eq. ()20|) is a linear differential equation in n, with a non- homogenous term $, depending 
on a parameter p. The general solution of the homogenous equation is: 

r(^)(«:;p) = r(^)(/co;p)e^-''o^^(^'), (22) 

where 

7(K) = (iV + 2)/f (k)^7(«:). 
14 



Note that, as k — 0, 7(k) ~ ^ with ^ > 0. The value of kq is to be chosen large enough for 
the LPA solution to remain a good approximation at this value of k (for example, kq = 10 p). 
A particular solution of the total equation can be found in the form: 

r(4)(/€;p) = fW(«:;p) (23) 

One gets 

T^'\k;p) = r^^K",p) e-^«o"^^("') + r(^)('^o;p), (24) 
so that the general solution of equation (j20|l can be written as 



r(^)(K;p) 



e^'^'^o^^^'^'). (25) 



This expression has the expected behavior: the last exponential in the r.h.s. guarantees that 
T^^\k,p) ^ as a power law when k 0. 

The behavior of the solution is shown in fig. |3 above and compared with the other 
expressions that we had for T^^\ for the particular values p = u and p = m/10. One can 
appreciate the effect of the improved procedure: the corresponding curve follows the direct 
LO one down to k ~ p/10, then it correctly goes to zero when k — > while the direct LO 
one does not. This analysis confirms the importance of keeping the feedback of the flow in 
the s-channel, that is, of solving the flow equation in this channel, in order to get the correct 
behavior at k ^ 0. 

As suggested by fig. |3J the range of values of k where r(^^(K;p, — p, 0, 0) obtained in the 
direct procedure exhibits its pathological behavior decreases as p decreases. In fact, as 
p — > 0, V^'^\k]p^ — p, 0, 0) ^ as K ^ (since the 4-point function is then given by the 
LPA'). One can argue at this point that the unphysical behavior of V^^^k) at small k has a 
moderate influence on the calculation of the NLO result for the self-energy. Indeed, as we 
shall see, it is mostly the region k ~ p that contributes significantly to the flow of TjNLoip)'-, 
in this region, the estimates of the 4-point function obtained in the direct and improved LO 
are almost identical (the difference with the initial ansatz is also small). This is why we have 
used the direct procedure in the first estimate of T^nlo presented in However, as we 
shall see in the next section, the improved procedure turns out to be much more accurate. 
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C. Improved calculation of T^^\K;p,—p,q,—q) 



The procedure described in the previous subsection only apphes to the 4-point function 
at g = 0. It is only for this value of q that we can write eq. as a closed equation: if 
g 7^ 0, the 4-point function in the l.h.s. and those in the r.h.s. are evaluated in different 
momentum configurations. However, we note that, in the calculation of T,nlo{i^]P), we need 
T^^\k]p, —p,q, —q) only for q < k, i.e., in a range of values of q that vanishes as k ^ 0. 
In that range, we expect T^^\k;p, —p, q, —q) to differ very little from T^^\k;p, —p, 0,0). 
Proceeding as in the previous subsection, we single out the s-channel and rewrite eq. as 

Kd^T^^'liK^p, -p, q, -q) = lf\K)V%{K;p, -p, 0, 0)rg(/€; g, -g, 0, 0) + $i2(«:,P, g) (26) 

where, now 

*i2(/t,P, q) = jf\K]p + g)f SJ^.(k;p, g, 0, -p - g)f ^^^.(k; -p, -g, 0,p + g) 



+ Jd (f^;?- (iWijif^; p, -q, o, -p + g)ry^.(K; -p, q,o,p- g) 

2^d \'^)'^12ll 

6i2^{K;p,q). 



- lld\^PSammi^;P, -P, 9, -9, 0, 0) (27) 



In eq. (EHj) . we use the fact that g < k to replace rlf-j{K; g, — g, 0, 0) by the initial ansatz (we 
have seen earlier that this is a good approximation): 

r;g.(«:; g, -g, 0, 0) ^ (iV + 2)6,,g{^). (28) 

As a result, all the g-dependence is now in the function $(k; p, g). For the other factor in the 
r.h.s. of eq. T\f-j{K;p, —p, 0, 0), we use the improved solution obtained in the previous 
subsection, which we shall denote as T^^2ii{'^'i Pj ~P^^}^) = Si2^'^{i^,p)- The 4-point function 
^uiii'^'^ p, —p, g, — g) can then be obtained by simply integrating the r.h.s. of eq. (j26|) between 
A and k and adding the initial value (see eq. (12)). One then gets 

rS,(«:;p,-p,g,-g) = 5i2(A^ + 2)^ (29) 



where 

T\k;p) 



r'^ riff' 

+ 5i2 / —[l{K')T\K';p) + ^K';p,q)], 

{N + 2)g{K) ifK > Ko 

r ^$(«:'; p) e/«' l^^'^C'") + {N + 2)g{^o) e^^^o "^'^^^'^ if < • 



(30) 
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FIG. 5: The function T^'^\k,,p, — p, 0,0) (in units of A) for the values p = lOu and p = u/10 as a 
function of k/u, for a = 0.6 and a = 0.9. The a-dependence is largest when p Kc — u/10. When 
p = lOu ^ Kc, r*^'''' is independent of a. 

D. The dependence on a 



In paper I we discussed the dependence of S^o on the parameter a that we introduced 
in one of the approximations (approximation A2) used to construct the initial ansatz for 
the 4-point function (see eq. (jllj) ). Recall that this approximation consists in replacing 
propagators such as G{p + q) in the r.h.s. of the flow equations by G'lpa'(q')0('^^ — a^p^). 
We found in paper I that Sio(p; «) obeys an approximate scaling law, T,Loij>\ — ^(ap). 
Here, we discuss the a-dependence of the LO expression of the 4-point function. For simplic- 
ity, we shall discuss the dependence of the expression obtained using the direct procedure 
(sect. Iir^ . The a-dependence of that obtained with the improved procedure is essentially 
identical. 

To study the a-dependence of the 4-point function, it is convenient to separate 
^mA'^'^P^^'P^^^"^) i^to three contributions: that of the s-channel (see fig. [T]), denoted 
by F*^^'"''!; the sum of the contributions of the t and w-channels (see fig. 12)), denoted by 
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p(4)[t+'u]. gj^ally the contribution of the 6-point function (see fig. El; denoted as r^^^^t^l. That 
is, r*^^^[*] denotes the contribution to the 4-point function obtained when only the channel i 
is included in the calculation of T^^^ according to eq. Q- There is an important difference 
between r(^)W and r^^)!^] 

on one side, and r^^^'*"*""] on the other side: while p flows through 
the loop in r*^^-"*'^'"] (see fig. 12)), in the other two cases it does not, so that the p-dependence 
of r'-'^-'t'*' and r*^^-*t^l comes entirely from the vertices. The latter are those of the initial ansatz 
r^^) and T'^^\ which depend on p and a approximately only through the product pa (see 
paper I, sect. Ill and app. E}. On the other hand, the p-dependence of r*^^)[*+"l is mainly 
due to the explicit p-dependence of the loop in fig. |21 which is given by the (a-independent) 
function jjf\K;p): Since J^p{K]p) is only important when p (see fig. 9 in paper I), the 
contribution of the t + u channels is non zero only in a region where the vertices in fig. |21 
are essentially the {p and a independent) LPA' ones (see e.g. fig. EJ. Then, one expects 
^mf~^^\'^''P^ ~Pi ^' t)e almost independent of a. 

Fig.Elshows that the total 4-point function V^^^ (k, p, —p, 0, 0) is in fact almost independent 
of a when p = 10m, which indicates that the contributions of the t and u channels dominate 
for this value of p. The same holds for all values of p much larger of much smaller than u. 
Only in the intermediate momentum region p Kc ^ u/10, is the variation with a important, 
which reflects the fact that in this intermediate range of momenta, the contributions of F*^^^'*! 
and r'^^^t^l are of the same order of magnitude as that of r*^^)[*"^"l, as we shall verify later (see 
e.g. fig. 7 below, and the related discussion concerning the a-dependence of SjVLo(p))- 

III. THE SELF-ENERGY AND A{ip^) AT NLO 

We have now all the ingredients to calculate the self-energy at next-to-leading order. 
Recall that the physical self-energy at criticality is given by (see eq. (1.108)) 

<5i2S(p) = ^jyK' J -0-^G\K';q)d,,R,iq) (T^Sui'^';p, -p,q, -g) - rgU'^'; 0, 0, g, -g)) . 

(31) 

In order to get J^nlo, one needs to insert in the r.h.s. of eq. (p?T|) the LO expression for the 
4-point function F^giK^SP' which has been calculated in the previous section, and 

the LO propagator given by: 

G-\k- q)=q^ + Sio(«:; q) + (32) 
18 
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FIG. 6: Self-energy at NLO (in units of A^) as a function of p/u for various values of a: a = 0.6 
(full line), a = 0.75 (dashed line) and a = 0.9 (dot-dashed line). 

with Slo(^; i) the LO expression of the self-energy, given by eq. (I. 111). 
A. Self-energy T^nlo 

In this subsection, we present numerical results obtained for = 3 and = 2 in order 
to illustrate the main features of the self-energy. We shall use here the LO estimate of the 
4-point function derived in the direct procedure. In the next subsection, we shall discuss 
results obtained with the improved procedure; we shall also present results for N ^ 2. 

Fig. El displays the self-energy at NLO for various values of a. It is to be compared with 
fig. 12 in paper I, for the LO behavior of S(p). Note, however, that in paper I S^o is plotted 
as a function of {ap/u) in order to exhibit its approximate scaling property T,lo{ol'-,p) ~ 
S(q;p). Here Satlo is plotted as a function of (p/u); as can be seen on fig. El ^nlo depends 
much less on a than S^o- In fact, both the IR and the UV regimes are nearly independent 
of a. It is only in the intermediate momentum range [p ~ Kc) that T^j^lq exhibits some 
dependence on a. 

At high momentum, one expects S(p) to be given by perturbation theory, that is, one 
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expects E(p) ~ \n(p/u). Recall that in paper I, we found that T,lo{p) indeed behaves in this 
way, but the coefficient in front of the logarithm differed from that of perturbation theory 
by 7%. As we have discussed in paper I, perturbation theory is recovered exactly at high 
momenta when one performs iterations. In particular, the 2-loop result is exactly reproduced 
at NLO. We have verified that the coefficient in front of the logarithm in T^j^fioip) large 
p, is correctly obtained, to within the numerical accuracy with which it can be determined 
(about 0.5%). 

In the IR regime, the power law behavior already reproduced in LO, p"^ + T,Loip) ~ 
is essentially not modified: however the numerical calculation is more involved in NLO, 
leading to a loss of accuracy that prevents us to determine the value of rj with any useful 
precision. 

As was the case in LO, the a-dependence of T^nlo is intimately connected with its mo- 
mentum dependence. Following the analysis that we did in sub-section III Dl to understand 
the variation of F'-'^-* with a, we split Tj^lo into three separate contributions: we define 
S'*], S^*"*""! and as the contributions obtained, respectively, when only F'^^^'^'l, F'^^^t*"'""! 
and F(4)[61 are included in the r.h.s. of eq. (j!?T| . The properties of the LO 4-point function 
discussed in the previous section imply that T}-^\a]p) ~ T!^^\ap), T}-^\a]p) ~ Y}-^\ap) (i.e., 
the same dependence on a as S^o (see paper I, sect. IV A)), while E'*"^"! is expected to be 
roughly independent of a. These properties are well verified in our numerical calculations. 

The three contributions of the self-energy are shown in fig. [7| together with their sum 
Satlo- While is positive, the two other contributions are negative. The latter property 

can be understood as follows: In calculating, say, S^*', one puts in eq. (pT|l only F'^^'^t*' 
which, in turn, is calculated with only the ffist term of eq. Q. In the latter, the flow 
is evaluated with the initial ansatz F^^^ which, as can be seen in sect. Ill C of paper I, 
verifies, in all regions of momenta, V^^\k]p^ —p^Q, -q) > f(^)(K;0,0,g, -g) > 0. Since the 
integration over k' that is needed to obtain F(^)[*l(/t) from eq. Q runs from A to k, we have 
F'^'')W(/t;p, — p, g, — g) < F(^)W(fi;; 0, 0, g, — g) so that the integrand in eq. (p?T|l is negative, 
yielding eventually S''^^ < 0. A similar analysis can be done for S^^^. Fig. [3 shows that, as 
was the case for the 4-point function discussed in the previous section, the t and u channels 
dominate, except in the intermediate momentum region {p ~ k^) where the contributions of 
the three channels are of the same order of magnitude. This, together with the a-dependence 
of St''!, E'^I and Sl*"*""] recalled in the previous paragraph, explains the behavior seen in fig.lHl 
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FIG. 7: Self-energy T^isfLoip) (™ units of A^) as a function oi p/u (triangle) at a = 0.83 and its 
three partial contributions (diamond), — Sl'*! (circle) and — S^^] (square). 

Fig. [7| also reveals an interesting feature of the present approximation, for which we have 
no simple explanation: to within numerical accuracy, T,^^^ and S^^J appear indistinguishable. 
This property remains true for other values of A^. 

B. Calculation of A{ip'^) 

We turn now to the calculation of the changes of the fluctuations of the field caused by 
the interactions: 

N ~ j (27r)3 + p2 J - 27r2 y p 

As recalled in the introduction, and more thoroughly in paper 1, this quantity is very sensitive 
to the intermediate momentum region and constitutes a stringent test of the calculation. 
In the following, we shall refer to the quantity in square brakets in eq. (jSBI), multiplied by 
l/(27r2), as to the integrand. Note that in the range of momenta where this integrand is 
significant (see e.g. fig.|Hlbelow), S(p) ^ p'^, so that the integrand can be well approximated 
by -S(p)/(27r2p). 



p2 + i;(p) 



P 



(33) 
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FIG. 8: The three curves represent the integrand of eq. (|33j) (in units of A) calculated with 
only (diamonds), Sl'^l (circle) and S^^l (squares) contributions to the self-energy at NLO, 

respectively, as a function oip/u (the points shown are those needed in the numerical calculation 
of the integral in eq. (|33j) ). The plots correspond to a = 0.83 and N = 2. 



The results for A(0^) will be discussed in terms of the parameter 



(C(3/2))^/^ Nu 



(34) 



The shift, caused by weak interactions, in the temperature of the Bose- Einstein transition 
of a dilute eras is 

As we have seen in fig. IHl T,nlo{p) is independent of the parameter a both at high and 
low momenta. However, in the crossover region {p ~ Hc) which determines the value of 
c, StvloIp) still depends on the value of a. It follows that the value of the coefficient c 
calculated with T,nlo{p) still depends on a. To understand better the a-dependence of the 
NLO predictions, one can write the coefficient c as the sum of three expressions, each of 
them containing only one of the three parts of the self-energy. The three contributions to 
the integrand yielding c are displayed together in fig. |H1 Because of the approximate scaling 
discussed above, both St*' and S^^l contribute to c with a term proportional to a, while the 
contribution of S'*^"] is essentially independent of a. This afine behavior of the NLO result 
for c is indeed observed in fig. IHl below. The negative slope is due to the negative sign of 



22 




0,6 0,65 0,7 0,75 0,8 0,85 0,9 

a 

FIG. 9: The coefficient c as a function of a for N = 2: LO (circles) from paper I, NLO direct 
(squares), NLO improved (diamonds). 

SWW and SW[61. 

The a-dependence remains a source of uncertainty in the calculation of c. As can be seen 
on fig. El when we move from the LO calculation to the direct NLO, to the improved NLO, 
the dependence on a decreases, and so does the corresponding uncertainty in the calculated 
value of c. We regard the variation in the value of c when a runs form .6 to .9 as a large 
overestimate of the uncertainty related to the choice of a. In fact we can eliminate much 
of this uncertainty by following a procedure suggested by the results plotted in fig. El since 
the curves representing c as a function of a have opposite slopes at LO and NLO, one can 
invoke a principle of fastest apparent convergence to choose as best estimate that given by 
the value of a for which the two curves cross: At this point indeed, the NLO correction 
vanishes. One thus obtains the value c = 1.44 (the crossing point being at a = 0.83) when 
we use the direct LO expression of the 4-point function, and c = 1.30 with the improved 
LO (corresponding to a = 0.77). The improved NLO calculation is thus in remarkable 
agreement with the lattice data: 1.32 ± 0.02 Q and 1.29 ± 0.05 
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c 


= 1 


iV = 2 


N = 3 


iV = 4 


N = 10 


iV = 40 


N = oo 


lattice [161 




1.32 ±0.02 












lattice 117} 




1.29 ±0.05 












lattice [181 


1.09 ±0.09 






1.60 ±0.10 








7-loops [19j 


1.07 ±0.10 


1.27 ±0.10 


1.43 ±0.11 


1.54 ±0.11 








large iV [22j 














c = 2.33 


this work 


1.11 


1.30 


1.45 


1.57 


1.91 


2.12 





TABLE I: Summary of available results for the coefficient c. The last line contains the results 
obtained in this work by using the improved LO approximation for the 4-point function. 

We have also repeated our calculation for other values of for which results have been 
obtai^ed With othe. techn^ues. e,th. the .att.ce technique UW, or va.atio„ally 

improved 7-loops perturtbative calculations These results are summarized in table I 
and fig. ^1 (other optimized perturbative calculations have also been recently performed, 



and are in agreement with those quoted here; see [20l|2l|). For small values (A^ ^ 10), our 
results fulfill all the numerical tests that we have described in this paper. For < 4, where 
we can compare with other results, the values of c obtained with the present improved 
NLO calculation are in excellent agreement with those obtained from lattice and 7-loops 



calculations 
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What happens at large values of A^ deserves a special discussion. As seen in fig. ^] the 
curve showing the improved leading order results extrapolates when A^ — oo to a value that 
is about 4% below the known exact result j^]- A direct calculation at very large values of 
A^ is difficult in the present approach for numerical reasons: since the coefficient c represents 
in effect an order 1/A^ correction (see it is necessary to insure the cancellation of the 

large, order A^, contributions to the self-energy, in order to extract the value of c. This is 
numerically demanding when A^ ^ 100. Fig. also reveals an intriguing feature: there 
seems to be no natural way to reconcile the present results, and for this matter the results 
from lattice calculations or 7-loop calculations, with the calculation of the 1/A^ correction 
presented in ref. |i23|: the dependence in 1/A^ of our results, be they obtained from the direct 
LO or the improved LO, appear to be incompatible with the slope predicted by the 1/A^ 
expansion. 
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FIG. 10: The coefficient c (obtained with the fastest apparent convergence procedure) as a function 
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of 1/A^. Our NLO results, are compared with results obtained with other methods: lattice 

ll^ (diamonds) and 7-loops perturbation theory (squares), all of them with their corresponding 
error bars, together with the iV — > oo result |22] (circle), and the extrapolation following the \/N 
correction calculated in ref. 

IV. SUMMARY AND OUTLOOK 

The quality of the results that we have obtained for the parameter c that characterizes 
the shift in the transition temperature of the weakly repulsive Bose gas is encouraging. It 
demonstrates that the method that we have developed in order to solve the Non Perturbative 
Renormalisation Group is capable indeed to yield the full momentum dependence of the n- 
point functions, in physical regimes where other approximation schemes are limited. 

It is of course difficult to quantify the size of the theoretical uncertainties in this approach. 
We have commented in the previous section about the source of uncertainty related to the 
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choice of the parameter a. A better estimate of the accuracy of the whole approximation 
scheme would be to perform one more iteration. However the numerical effort involved in 
the calculation of the relevant multidimensional integrals is non negligible, and furthermore 
the approximate equations that need to be solved in order to get the initial ansatz for the 
n-point functions become increasingly complicated (see for instance app.El- Thus it seems 
unrealistic to imagine doing a further iteration, going say to the next-to-next-to-leading 
order. What remains then, as a measure of the quality of the approximation scheme, is the 
direct comparison with results obtained by other, reliable, non perturbative methods, such 
as lattice calculations. 

We have however started exploring an alternative approach which may have, among other 
features, the capability of yielding estimates for theoretical uncertainties. This approach 
builds on the present approximation scheme, but brings to it conceptual simplification. 
As we have seen when discussing results, both in paper I and in the present paper, the 
approximation A\ introduced in paper I is accurate in most situations encountered when 
solving the flow equations. This approximation assumes that the vertices in the r.h.s. of 
the flow equation are smooth functions of the momenta and exploits the fact that the loop 
momentum is bounded to neglect some of the momentum dependence. Once this is done, 

n 

as shown in p^, one can relate directly the higher n-point functions that arise in the r.h.s. 
of the flow equations to derivatives of the ra-point function whose flow is being studied with 
respect to a constant background field. One then obtain closed equations. This method 
allows us to bypass both the approximation ^3 needed for n-point functions of high order, 
and also the approximation A2 used to implement the decoupling of high momenta in the 
propagators. The price to pay is that the new equations are not only differential equations 
in the variable k, but also partial differential equations with respect to the background 
field. However, they can be solved numerically, as will be demonstrated in a forthcoming 
pubhcation jlj]. 

APPENDIX A: THE INITIAL ANSATZ FOR rf2l,„„(K;p, -^j, 0, 0, 0, 0) 

In order to obtain the explicit expression of the initial ansatz for F^^-* we follow the same 
steps as in paper I, sect. HI, when constructing the initial ansatz for F^^-': We use the three 
approximations A2 and ^3 to get an approximate equation for the flow of T^^\ which 
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we then solve semi-analytically. 

The exact flow equation for T^^^ can be written in the following form, exhibiting three 
kinds of contributions: 

i^d^iL^e = TTG\d,R^[-r%Grfl,Grl^l^ + UpeTmntations 

+ ^illsAjGT^i + 14 permutations - ^rfaLseu | • (Al) 

We start by implementing approximations Ai and A2- To keep the discussion simple, 
we do so explicitly only for some typical terms of eq. (jAlj) . Take for example the following 
contribution involving the product of three T^^^: 

-TTG'nd^R^T%GT%,GT^^l^ = - j ^,ndMq')G{q')T%{q,PuP2, -q - Pi - P2) 
X G{{q + pi+ P2)^)rjj4fc(g + Pi+ P2,P3, Pa, -q + P5+ Pg) 

X ^((g-ps -P6)^)r26i(g-P5 -P6,P5,P6,-g)- (A2) 

After setting the external momenta to their values (ps = P4 = Ps = Pe = 0,p2 = — Pi), 
imposing g = in the vertices (approximation ^1) and replacing G{q + p) by Glpa', 6(1 — 
a^p^ /q^) (approximation A2), one gets: 

Trr2 o p p(4) ^p(4) ^p(4) _ 

0' o)r3S.(o, 0, 0, o)r(l(o, o, o, o), (A3) 

where we have also made use of the symmetry of the bosonic n-point functions. A similar 
contribution corresponding to a different permutation reads: 

-I^f\n)Qil - ay/K')T%{p, -p, 0, 0)T%,ip, -p, 0, 0)1^^(0, 0, 0, 0). (A4) 

All 45 permutations containing three T^^^ reduce either to the forms ()A3|) or ()A4|) : those 
where both the external legs carrying the non-zero momenta {p and —p) belong to the same 
r*^^) are of the type ()A3|) . whereas those where the two legs with non-zero momenta belong 
to two different T^^^ are of the type ()A4|) . 

The second kind of contributions is that which involve one F'^'^^ and one r(6). Among 
the 15 contributions of this kind, we have three different cases, depending on how the two 
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external legs carrying non zero momenta are attached: both on T^^\ both on T^^\ one on 
r*^^) and the other on T^^\ We explicitly write here one contribution of the latter type: 

Tr G\d.R.Tfl,^^GT^^l, = /f (/s:)e(l - aY / K')Tf^,,,,{p, -p, 0, 0, 0, 0)T%M 'P. 0, 0). 

(A5) 

Finally, the third type of contributions is that which involves F*^^-': 

- ^ Tr G^Kd,Rj:fi,,,,^^ = -^/f («:)Ff2L56..(P, 'P, 0, 0, 0, 0, 0, 0). (A6) 

Observe that while all expressions in eqs. ()A3|) and ()A4|) are known (the explicit form of 
F*^^^ can be found in paper I, sect. IllC), in the r.h.s. of eq. ()A5|1 appears the function F'^^^ the 
variable of the differential equation (jAlH . Evaluating all the 15 permutations which includes 
this function T^^\ one verifies that it appears either in the form T^^\p, —p, 0,0,0,0), or in 
the form F(^^(0, 0, 0, 0, 0, 0). The latter being simply the (known) LPA' expression (see paper 
I, sect. IIC), one ends up with a differential equation for the function T^^\k; p, —p, 0, 0, 0, 0). 
In order to solve it, one needs an initial ansatz for F^^^ that appears in ()A6|) . 

To get the latter, we follow approximation ^3. Let us write the LPA' equation cor- 
responding to ()A1|) . This can be obtained by deriving three times with respect to p the 
equation for the effective potential, and then setting p = 0. One gets: 

d^V'"ip = 0) = J -0^^dMQ)G'iq){-3iN + 2Q)G\q)iV"ip = O)f 



+ 3{N + U)Giq)V"'{p = 0)V"ip = 0) - ^^^V^W(p = 0) 



(A7) 



Defining (see also paper I, sect. 11 C) 

V^'^ {p = 0)=L = Kfzy-^'L (A8) 

we transform eq. ()A7|1 into: 



Kd^K = -m + 2Q)Id\^)9l + 3(Ar + U)gjji'\t,)K + 6)/f («:)/«. (A9) 

On the other hand, as suggested by the large limit (see paper I, sect. II D), one expects 
that once approximation (^2) is performed, the term containing F*^*^) to be proportional 
to the other ones, the coefficient depending only on k. However, while in the case of the 
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FIG. 11: The function F' function of n/u. 



equation for F*^^^ the three contributions involving products of F^^^ in eq. where identical 
at zero external momenta (thus giving a unique proportional factor F^), here we have two 
different contributions (those with three F*^^) and those with one F*^^) and one F^^-*). Thus, 
the way the contribution of F^^) can be distributed over the other two is not unique. To 
remove the ambiguity, we distribute the various terms as they appear in the LPA', eq. ()A9j) : 

1 



(iV + 6)/f (/s:)/. = 3(iV + 26)P;> {K)F',gl - 3{N + U)I^'> {K)F^g^K, (AlO) 



r(4)/ 



(3), 



where F^ is the function defined in eq. (1.44), while eq. ()A10|) can be taken as the definition 
of F^. Using eqs. (1.42), (1.44) and (jX8|) . one then sets: 



F' 



2\2 



{N + IA){N + A)hl N + Q 



2{N + 



(All) 



The function F'^ is shown in fig. ^2 

We are now in the position to perform the approximation ^3 in the flow equation of V^^\ 
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i.e., in eq. ()A1|) . This amounts to the replacement: 



- - Tr ^ rf^,,,,sG'd^R^] ^Tr {c'd^R^ {K^%Gr%,Gr^l, + 44 permutations 



3{N + U)F^r[f2^^jGrf^>^^ + 14 permutations 



(4) 



(A12) 



At this stage, we have all the ingredients to write, and solve, the (approximate) equation 
for r^^''. After rewriting eq. ()A1|) with the use of eq. ()A12|) . and evaluating all the terms as 
in eqs. ()A3|) . ()A4|) and ()A5|) . one ends up with an ordinary differential equation where the 
dependent variable is T^^\k;p, —p, 0,0, 0,0). To write explicitly this equation it is useful 
to use the fact that rf2345e{p, —p,0,0,0,0) is completely symmetric under the exchange of 
indices 3, 4, 5 and 6. Then, one can make the decomposition: 

r 123456 

(p, -p, 0, 0, 0, 0) = 1(6)512(534556 + 535^46 + 536^45) 
+ rf ) [524(513^56 + 5i5536 + 5i6535) + 5i4 (523556 + 525536 + 526535) 

+ 5i3525546 + 5i3526545 + 5i5523546 + 5i5526543 + 5i6523545 + 5i6525534)] , (A13) 

Finally, taking the trace over the tensor indices, and doing a lengthy, but straightforward 
calculation, one gets: 

«^5.ri6) = (1 - F^)if\K) U \i{N + 8)rf + 4rf ) 



N + 6 + 



N+2 2 



NJ \gap 



- 3(1 - F'^)l!^\^)glg^, i + 10 + ^) { ^ 



N \gap 



(A14) 



(6) 



;i - F,)lf\K) { g,{N + 16)r^^^ + 2Map ( — 

^9oip 



.(6) 



9k 



2 

JV+8 



- m-F'^)lf\n)glg^,[^ 

9ap 



9k 



2 

JV+8 



(A15) 



when ap > k. When ap < k, one can show that the LPA' expression for the 6-point vertex: 

Tfl^, = 5,2KiN + 2)iN + A) iap<n), (A16) 



30 



is a solution of the equation that one gets. Then, as it is a first order differential equation 
in K, the expression ()A16|) is the solution for the case ap < k. 

Returning to the case ap > k, one can diagonalize eqs. ()A14|) and ()A15|) . the eigenvalues 
and eigenvectors being 



Xa = 2{N + 8) 
Xb = N+16 



One thus can write: 




with eigenvector (1, 0) 
with eigenvector {—A/N, 1) 



-A/N 




(A17) 



(A18) 



where verifies the equation: 

nd^aa = (1 - (/€)2(iV + 8)g^aa 

N+2 

+ (1 - F^/f («:) ( AT + 6 + 1^ ) Kg^, (^^ ' 



-3(1 - ( AT + 10 + ^) {g^,r 



, N+6 
' N+8 



(A19) 



while the equation for aj, is simply the equation for F^^-* ()A15|) . 
Both equations ()A15jl and ()A19|1 have the following form 



its solution being 



A 

T{n)=T{ap)[^]^^\ [^K^f^ 

^9apJ Jap \9k,' 



9k 



A 

JV+8 



(A20) 



(A21) 



Unfortunately we could not succeed to solve analytically the equation above, as we did 
for the initial ansatz for F*^^^; there, the key equation was (1.43), but we could not find a 
similar one here. 

Imposing the continuity condition, in k = ap, dictated by eq. (jA16j) : 

T^^\ap) = T^^\ap) = K 



''ap 



(A22) 



which gives 



aa{ap) = 
ab{ap) = hap 



l + J^]hap 



(A23) 
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the solutions for aa{K,) and ah{n) are 

2 



log op/ A 



iV+2 
iV+8 



3(1 - F:,)I^'\^){N + 10 + ^) 

ab{K) = hap 



oJV+6 



iV+16 



9ap 



N+ia f , 2 



log op/A \9kJ I \9apy 



- 12(l-F^)tf)(«:)^?^^7.p(f^ 



2 

JV+8 



(A24) 



Finally, when k < ap, one has 

rf4,(P, -P, 0, 0, 0, 0) = 6uN{N + 2)aa{t^) (/t < ap). (A25) 

APPENDIX B: PRODUCTS OF FUNCTIONS f^^) IN EQ. © 
1. The s and i-channel contributions 

In this appendix we obtain explicit expressions for the products of functions F^^^ that 
appear in the r.h.s. of eq. Q, where F*^'^-' is the initial ansatz for the 4-point function. 
All the needed expressions for the 4-point fucntions that are needed here are those for 
T^^\pi,P2, 0, —pi — P2) that can be found in paper I, sect. Ill B. 

1) The s-channel contribution 

Here we consider the product T^iji'^'iPi 0) ^)^uiji'^'i 'ii 0) 0)) where, because of the 
regulator, q < k. There are two regions to examine: 
a) ap < K, 

Both 4-point functions are in the region (a) of paper I, sect. Ill C. After a simple 
calculation one gets: 

f SS,(P, -P, 0, 0)f!£.(g, -g, 0, 0) = {N + 2fgl5,2. (Bl) 
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b) ap > K 

Here, while t\fliq, -q, 0, 0) is still in region (a), the other vertex is in region (b) of paper 
I, sect. Ill C. In this case one gets: 

JV+2 

f SS,(P, -P, 0, 0)f[g(g, -q, 0, 0) = (iV + 2fg,g^p (^) 5u. (B2) 
2) The i-channel contribution 

Here we consider the product f —p, 0, 0)f[f^j{q, —q, 0, 0). As q < k, and a < 1 (and 
thus aq < k), one has two cases to study: p > \p + q\ and |p + q'| > p- 
A) p>\p + q\ 

In paper I, sect. Ill B, we assumed Pi > P2 > \pi + P2\- Using the symmetry of the 
bosonic n-point functions, we can conveniently rewrite the product as 

^uiM ^' °' ~P ~ ^)^2iiji-P^ O'P + = ^ijiiiP^ -P - 1^ 0' 'l)^^jiiiP^ -p - o> 

(B3) 

and the expressions of the vertices to consider are those of either regions (a), (b), or (c) of 
paper I, sect. HI B, depending on the value of k (region (d) never enters, because q < k): 
a) K > ap 

The two vertices are in region (a) of paper I, sect. Ill B. The product is simply: 

f J -Up, -P - q, 0, g)f g;,(p, -p - q, 0, q) = 3{N + 2)g%,. (B4) 
h) ap> K> a\p + q\ 

Now, both vertices are in region (b) of paper I, sect. HI B. A lengthy but straigtforward 
calulation yields: 

TykP' -P - 0' (l)^fjii{P^ -P - 0, q) = 

a-^— o.^vra^)(i^)- 

+ (^^N' + 6iV + 30 - ^^^±l^VN^ + 4N + 20^ (j^^ ""^^ I 6u 

(B5) 

c) a\p + q\ > K 
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Both vertices are now in region (c) of paper I, sect. Ill B. After another straigtforward 
calulation one obtains: 

f hL(p> -p - 0' 9)rSzb, -p - q, o, g) = 



^1 . . ^2 

N+S / n \ 



iv I +b- ' 



9a\p+q\) \P+l\ Y 



a\p+q\ 



9k 



9a\p+q\ 



N+8 



^2 

N+8 / n \ 'v+g 



^+ ( 1 +6^ (9k 

\9a\p+q\/ \ya\p+q\ , 

N'^ / n \ o 



(B6) 



where &„|p+g| and r^i^^^i follow from eqs. (1.92) and (1.93) respectively. 

B) p < \p + q\ 

In this case, we reorder momenta and indices as: 

i'liijip, Q, 0, -p - q)i'2Hj{-p, -q, o,p + q) = i'fiui-p - q,p, o, q)i'f2ii{-p - q,p, o, g). (B7) 

Similarly as in previous case (A) one has three regions, k > (y\p + q\, (y\p + q\ > k > ap 
and ap > k. For each region the result is the same as those in eqs. ()B4|) . ()B5|) and ()B6|) . 
but exchanging p with |p + g|. 
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by essentially the LPA' vertices, and the solution of the LPA' is also very approximate. Fi- 
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